Precision spectroscopy of the hydrogen molecular ion HD + 
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Expectation values of the Breit operators and the Q terms are calculated for HD + with the 
vibrational number v — — 4 and the total angular momentum L — — 4. Relativistic and radiative 
corrections to some ro-vibrational transition frequencies are determined. Numerical uncertainty 
in RocCt 2 order correction is reduced to sub kHz or smaller. Our work provides an independent 
verification of Korobov's calculations [Phys. Rev. A 74, 052506 (2006); 77, 022509 (2008)]. 
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I. INTRODUCTION 



II. THEORY 



Hydrogen molecular ions, such as H^ and HD + , can 
be used [l|, [|| to deduce an improved value of the proton- 
electron mass ratio by comparing experimental and the- 
oretical spectroscopic data. To this end, several experi- 
ments have been setup for measuring high precision tran- 
sition frequencies in HD+ ||, [| and H+ [|. For HD+ 
in particular, the (v,L) — (4,3) — > (0,2) transition fre- 
quency has been measured at the 2 ppb level Q . In order 
to reduce the current uncertainty of 4.1 x 10~ 10 @ in the 
proton-electron mass ratio, both experiment and theory 
should reach a precision of sub-kHz or better. 

For light systems, nonrelativistic QED (NRQED) @,H 
approach is used to expand the energy in powers of the 
fine structure constant a. Nonrelativistic energies of H^ 
and HD + have been variationally calculated at the preci- 
sion of 10~ 15 @, IH for a wide range of vibrational states 
and at the precision of 10 -30 [10b14| ] for some low vibra- 
tional states. Relativistic and radiative corrections of or- 
der Rood 2 , i?ooa 3 , RooCt 4 , and RooCt 5 have been system- 
atically calculated [H|,[l6| for Hj and HD + ro-vibrational 
states (v = — 4, L = — 4). The leading order R^a 2 
relativistic corrections are now available with very high 
precision for the H 2 vibrational states (0,0), (1,0), and 
(0, 1) 17] . Recently, improved values of R^a 3 order cor- 



rections have been achieved by_ 
of the Bethe logarithm term |l£ 



an intensive calculation 



The purpose of this Brief Report is to present our in- 
dependent calculations of the Breit operators and the Q 
terms for the HD + ro-vibrational states (v = — 4, L = 
— 4) . Some important ro-vibrational transition fre- 
quencies are determined, which provides a verification 
of Korobov's theoretical work. In addition, the numer- 
ical results presented here might serve as a benchmark 
for other theoretical methods. We use atomic units 
(h = e = m e = 1), unless otherwise stated. The fun- 
damental physical constants involved are taken from the 
2010 CODATA recommended values d. 



Consider the hydrogen molecular ion HD + . After sepa- 
rating the center of mass coordinates for the system, the 
eigenvalue problem for the nonrelativistic Hamiltonian 
H becomes [l4[ 



H ^ = E a ^ , 



H 



A 12 V ri • V,. 2 + V , 



(1) 
(2) 



where i"i and r 2 represent respectively the position vec- 
tors of the electron and proton, relative to the deuteron 
situated at the origin, Ai = —(1 + m c i)/(2m c i), A 2 = 
-(l/m d + l/m p )/2, A12 = -l/m d , V = -1/n + l/r 2 - 
l/ri2 is the Coulomb interaction, and r 12 = i"i — r 2 . The 
energy eigenvalue problem for Hq is solved variationally 
using the basis set in Hylleraas coordinates 



fcifcfa.ra) = r\r^ 2 e- ar ^y^ (r u r 2 ) , (3) 

where yff^ (ri , r 2 ) is the vector coupled product of spher- 
ical harmonics for the electron and the proton to form 
a common eigenvector of L 2 and L z . More details on 
the construction of basis set for HD + may be found 
in [l2|, [3| . It should be pointed out that this basis set 
differs from the one used by Korobov [l5|, [TB| ■ The basic 
type of integrals required in the calculation of matrix el- 
ements can be evaluated analytically |19| using Perkins' 
expansion for r\ 2 . The procedure for handling singular 
integrals that appear in the evaluation of Breit operators 
can be found in [20(. 

The leading-order relativistic corrections due to the 
Breit operators are well established, which can be found 
in Refs. [2l| - [23| . The complete spin-independent part of 
order R^a 2 term is 



=a 2 (H Brcit ) 



(4) 
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In the above, the Darwin term 7r/(2rop)£(ri2) is the 
nuclear spin dependent recoil correction for the spin-^ 
particle, such as proton. This term vanishes in case 
of spin-0 or spin-1 nucleus, such as the 4 He nucleus or 
deuteron [24| . It is noted that <5(r2) is virtually zero due 
to the molecular nature of the system. 

Furthermore, the spin-independent radiative correc- 
tion of order R^a 3 may be expressed as (25l-[27|: 
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where j3(v,L) is the Bethe logarithm, Q{vi) and 
are the Q terms introduced by Araki and Sucher [2 




Q(r) 
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4 7rr 3 



(lnp + 7£ ;)(5(r; 



(7) 



and 7£ is the Euler gamma constant. 

For higher order corrections, such as orders RooCt 4 and 
RooCt 5 , we follow the work of Ref. [161 ] . Thus, RooO 4 order 
non-recoil relativistic and radiative corrections may be 
express as follows [To) : 



41 1 

as — 



2179 3523 



L 7T [ (348 804 2 4 J1 

x(5(r 1 ) + <5(r 12 ))+£; i ( c 4 1 ) |, (8) 

where -E^j is the Rooa 4 order relativistic correction. 

Since the electron is almost bounded in the ground 
state of hydrogen molecular ion, its wave function can 
be approximately expressed as a linear combination of 
two hydrogen- like wave functions ip e (*e) = C[ipia( T i) + 
ipisfoiz)]- Therefore, the most important Rood 5 order 
correction can be estimated using this approximate wave 
function EH H! 



1 



E^=a 5 -\n\ ., 
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where the constants ^6i, Aqq, and B50 are taken to be 
the constants of the Is state of the atomic hydrogen, 



TABLE I: Convergence study of the expectation values of 
<5(ri) and <5(ri 2 ) for the HD + (v = 4, L = 4) state, where N 
is number of basis functions. 



TV 


«(r0 


*(ria) 


2789 


0.188 729 941 487 


0.188 382 948 5531 


3415 


0.188 729 516 647 


0.188 382 898 6629 


4130 


0.188 729 943 378 


0.188 382 886 4800 


4940 


0.188 729 938 954 


0.188 382 884 0699 


5851 


0.188 729 935 696 


0.188 382 876 7837 


6859 


0.188 729 935 329 


0.188 382 876 7442 


8000 


0.188 729 935 332 


0.188 382 876 7368 


9250 


0.188 729 935 329 


0.188 382 876 7376 


Extrp. 


0.188 729 935 323(7) 


0.188 382 876 7375(1) 



-30.924- ■■ !3() . and 



i.e., A 61 = 5.419- •• [2 1 
B 50 = -21.556- ■■ jUJ. 

In addition to the relativistic and radiative corrections, 
one also needs to consider the contribution from the finite 
nuclear charge distribution. The leading-order correction 
is 



f -He 

^nuc " ^ 



ao J \ a o J 



, (10) 



where R p = 0.8775(51) fm and R d = 2.1424(21) fm [| 
are the root-mean-square charge radii of proton and 
deuteron respectively. 



III. RESULTS 

With the Hylleraas-type basis set of Eq. ((3]), the wave 
functions along with the corresponding nonrelativistic 
energies are obtained by solving Eq. ([1]) variationally. 
Then the expectation values of the Breit operators can be 
evaluated. In particular, the global operator method [32| 
is applied to the evaluations of V 4 2 , (V ri + V r2 ) 4 , 
<5(ri), and (5(ri 2 ). As an example, Table U shows a con- 
vergence study for the expectation values of <5(ri) and 
<5(ri2). One can see that an accuracy of about 11-12 sig- 
nificant figures is achieved for the most difficult state of 
(4, 4), where the nonrelativistic energy is calculated only 
to 16 digits. For ro- vibrational states (v — 0—4, L = 0—4), 
numerical results of the Breit operators are presented in 
Tables [U and Hill together with a comparison with Ko- 
robov's values [15|. Results for the Q terms are listed 
in Tables |V] and IIVI All expectation values are in good 
agreement with Korobov's values, although our results 
are more precise. The largest size of basis set used here 
is about 9000. 

We summarize the contributions up to Rooa 5 to two ro- 
vibrational transition frequencies in Table IVI1 where the 
values of the Bethe logarithm are taken from Ref. [l8| and 
the Rood 4 order relativistic correction taken from 
Ref. [la ]. For the (1,0) — > (0,0) transition, the numer- 
ical uncertainty in A.E a 2 has been reduced from 1 kHz 
in Korobov's value to the present 8 Hz, which is due en- 
tirely to the uncertainties in the fundamental constants. 
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TABLE II: Expectation values of 8(n), 8(r 12 ), , V 4 2 , and (V n + V r2 ) 4 for HD+ with v = - 4 and L = - 4. Korobov's 
results [L5J are listed in the second entry of each ro-vibrational state. 
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The correction AE a 3 to this transition has been obtained 
in Ref. [18| and reproduced here. The recoil correction 
of i?oo a 4 (m/M) and higher contributes at the level of 
relative 10~ 10 — 10~ n , which causes a theoretical uncer- 
tainty of 1 kHz in AE a 4. The largest uncertainty for 
the transition (1,0) — > (0,0) comes from the theoretical 



uncertainty of AE a s . The uncertainty in A.E nuc is due 
to the uncertainties in the proton and deuteron charge 
radii. For the transition (4, 3) — > (0, 2), both experimen- 
tal Q and theoretical [16| results are available. In our 
calculation for this transition, although the uncertainties 
in AE a -i and AE a 3 have been reduced to sub kHz, the 
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TABLE III: Expectation values of Rd p , Rde 



and R pe for HD + with v — 



4 and L = — 4. Korobov's results MM are listed 
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15.7580 


1.14792 


1 ~i a r\ A c\ 

1.14042 


(2,2) 


24.13150123361982(4) 


1.12583357622946(1) 


1.11536372709830(3) 




24.1315 


1.12583 


1.11536 


(3,2) 


31.275189712792(3) 


1.1054236322514(2) 


1.0923565391832(2) 




31.2752 


1.10542 


1.09236 


(4,2) 


37.275607599566(2) 


1.086630566009(1) 


1.07130988884988(7) 




37.2756 


1.08663 


1.07131 


(0,3) 


6.74277630617217739(6) 


1.16907791520748(2) 


1.1645223841766(2) 




6.74278 


1.16908 


1.16452 


(1,3) 


16.35924480624(3) 


1.145392790044(2) 


1.13750662824(1) 




16.3592 


1.14539 


1.13751 


(2,3) 


24.6527816700128(9) 


1.1234560568468(8) 


1.112634817630(2) 




24.6528 


1.12346 


1.11263 


(3,3) 


31.72221152212(6) 


1.103194064780(2) 


1.0898050864(1) 




31.7222 


1.10319 


1.08981 


(4,3) 


37.65328339(4) 


1.0845441798(5) 


1.068929626(1) 




OT /"'COO 

37.6533 


1 r\ O A f A 

1.08454 


1.06893 


(0,4) 


7.640242669140628(6) 


1.16553169532956(5) 


1.1604316913454(5) 




7.64024 


1 1 crro 

1. 16553 


1.16043 


fl 4) 


17. 1426952788899(3) 


1.1420581168039(4) 


1.133671151272(3) 




17.1427 


1.14206 


1.13367 


(2,4) 


25.3308170266259(7) 


1.120324362118(2) 


1.109043458742(1) 




25.3308 


1.120320 


1.10904 


(3,4) 


32.30232200847(1) 


1.100257957806(8) 


1.086448221081(5) 




32.3023 


1.10026 


1.08645 


(4,4) 


38.14194609(1) 


1.0817974074(2) 


1.0657990275(3) 




38.1419 


1.08180 


1.06580 



total uncertainty in the transition frequency remains as 
large as 70 kHz, which is from the AE a 5 term. 

In summary, we have presented an independent cal- 
culation of the Breit operators and the Q terms for the 
HD + ro-vibrational states (v — — 4, L = — 4), which 
provides a verification of previous theoretical results. 
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TABLE IV: Numerical values of Q(ri) for HD + with v = — 4 and L — — 4, where the second entry lists Korobov's results [l5| ]. 







v = 


v = 1 


v = 2 


w = 3 


v = 4 




n 
u 


-D 1 ^48fi997fif : ;fl81 QfnC9l 


-0 1 "391 1 "^^l 9Qf1 


-fl 1 90^1 087^71 971 ( 1 


-fl 1 971 fiQ9Q044 c i c il'1 


-fl 1 94Q c i74Q c i c i77^Q('^ , | 






-0.13486 


-0.13211 


-0.12955 


-0.12717 


-0.12496 


/ _ 

/ . — 


1 

1 


u. io4; f 4ZUU1UJ f oyz^jj 


-0 1 ^1 QQQfil 1 4'}Qfi7l'9 > l 


— D 1 9Q444 c ;978fi8 c ;f4^ 


— fl 1 97flfi78Q778n9('9^ 


— fl 1 948R1 Q897849l'9^ 






-0.13474 


-0.13200 


-0.12944 


-0.12707 


-0.12486 


L = 


2 


-0.13450284686085(8) 


-0.1317734666211(3) 


-0.1292308990657(3) 


-0.126866353373(3) 


-0.1246721506275(1) 






-0.13450 


-0.13177 


-0.12923 


-0.12687 


-0.12467 


L = 


3 


-0.1341475653900(1) 


-0.131437564457(6) 


-0.12891364214(1) 


-0.12656709998(4) 


-0.124390349(3) 






-0.13415 


-0.13144 


-0.12891 


-0.12657 


-0.12439 


L = 


4 


-0.1336801802328(6) 


-0.130995769843(4) 


-0.128496470584(6) 


-0.12617370653(3) 


-0.1240200070(4) 






-0.13368 


-0.13100 


-0.12850 


-0.12617 


-0.12402 



TABLE V: Numerical values of Q(ri2) for HD + with v — — 4 and L = — 4, where the second entry lists Korobov's results [TBI ]. 







v = 


v = l 


v = 2 


v = 3 


v = 4 


L 


= 


-0.1345911955685105(1) 


-0.131838977797(3) 


-0.129272929986259(7) 


-0.1268839283180(2) 


-0.12466388687278(6) 






-0.13459 


-0.13184 


-0.12927 


-0.12688 


-0.12466 


L 


= 1 


-0.13447097875759(1) 


-0.131725268104580(5) 


-0.1291654760689(3) 


-0.1267825085601(1) 


-0.1245683078029(8) 






-0.13447 


-0.13172 


-0.12916 


-0.12678 


-0.12457 


L 


= 2 


-0.13423194023120(3) 


-0.1314991894516(1) 


-0.1289518569532(7) 


-0.126580907975(1) 


-0.124378340936(7) 






-0.13423 


-0.13150 


-0.12895 


-0.12658 


-0.12438 


L 


= 3 


-0.133876827784(1) 


-0.13116338201(2) 


-0.12863460989(3) 


-0.1262815656(3) 


-0.12409633(1) 






-0.13388 


-0.13116 


-0.12863 


-0.12628 


-0.12410 


L 


= 4 


-0.133409659744(2) 


-0.13072170536(5) 


-0.12821744286(1) 


-0.1258880438(2) 


-0.12372570179(5) 






-0.13341 


-0.13072 


-0.12822 


-0.12589 


-0.12373 
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(4,3) -> (0,2) 
Present 


(1,0) H 


• (o, o) 


Present 


Korobov 116] 
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57 349 439.973 34(14) 


57 349 439.9717 


AE a , 


3 411.702 93(4) 


958.276 694(8) 


958.277(01)" 


AE a3 


-891.610 9(1) 


-242.126 26(4) 


-242.125(02) 


AE a , 


-6.457(1) 


-1.748(1) 


-1.748 


AE a , 


0.388(74) 


0.105(20) 


0.105(19) 


AE nuc 
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-0.125 2(17) 


-0.125(2)" 


AE tot 
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